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Abstract 

We provide a mathematical analysis of appearance of the concentrations (as Dirac masses) of the 
solution to a Fokker-Planck system with asymmetric potentials. This problem has been proposed 
as a model to describe motor proteins moving along molecular filaments. The components of the 
system describe the densities of the different conformations of the proteins. 

Our results are based on the study of a Hamilton- Jacobi equation arising, at the zero diffusion 
limit, after an exponential transformation change of the phase function that rises a Hamilton- 
Jacobi equation. We consider different classes of conformation transitions coefficients (bounded, 
unbounded and locally vanishing). 



Key words. Hamilton-Jacobi equations, molecular motors, Fokker-Planck equations 
AMS Class. Numbers. 35B25, 49L25, 92C05 

1 Introduction 

A striking feature of living cells is their ability to generate motion, as, for instance in muscle contrac- 
tion already investigated theoretically in the 50's ([H]). But even more elementary processes allow 
for intra-cellular material transport along various filaments that are part of the cytoskeleton. These 
are known as "motor proteins". For example, myosins move along actin filaments and kinesins and 
dyneins move along micro-tubules. In the early 90's, it became possible to device a new generation of 
experiments in vitro where both the filaments and the motor proteins are sufficiently purified. This 
lead to an improved biophysical understanding of the biomotor process (see, for instance, [T| [T5 | [23 | [TT] . 
and the tutorial book [T7|) and gave rise to a large cellular biology literature. The experimental obser- 
vations made possible to explain how chemical energy can be transformed into mechanical energy and 
to come up with mathematical models for molecular motors. The underlying principles are elementary 
and represent in fact the common basis for all biomotors. On the one hand, the filament provides for 
an asymmetric potential (this notion was introduced in the earliest theoretical descriptions by Huxley, 
[18)). sometimes referred to as the energy landscape. On the other hand, the protein can reach several 
different conformations. This can be ATP/ADP hydrolysis but five to six different states of the protein 
could be involved during muscular contraction. 
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In this paper we consider the following model: Molecules can reach / configurations with density, 
for each i = 1,2..., I, n^. A bath of such molecules is moving in an asymmetric potential seen diff'er- 
ently by the / configurations denoted, for i = hy ipi. Fuel consumption triggers a configuration 

change among the different states with rates Vij > 0, for i,j = 1,2..., I. Diffusion, denoted below by 
a, is taken into account. 



Thes simple considerations lead to the following system of elliptic equations for the densities 

("-i)l<i</- 

-cr-§^ni - ^(VV'i Ui) + vam = ^ VijUj in (0, 1), 

j^i (1) 
a^ni{x) + \/'ipi{x) ni{x) = for x = or 1. 

The zero flux boundary conditions means that the total number of molecules, in each molecular 
state, is preserved by transport (but not by configuration exchange). 
Throughout the paper we assume that, for i = 1, 

ni>0 in [0,1]. (2) 

The zero fiux boundary condition, motivated by the additional modeling assumption that total 
density is conserved, leads to the condition that, for all i = 1, 

Several biomotor models, including the one described above, were analyzed in [71 [U [191 [16] through 
optimal transportation methods. In [7J it is proved that there is a positive steady state solution that 
can, for instance, be normalized by 

/ y2 ^ii^)dx = 1. (4) 
■^0 1<^<I 

The simplest way to explain this fact is to observe that the adjoint system, 
-(y-^<Pi + ^ipi ■§^<Pi + i^ii<Pi = Y in (0, 1), 



(5) 



V dx 



^^--^ in {0,1}, 



admits the trivial solution (pi = (j)2 = ... = (pi = 1. This yields that is the first eigenvalue of the 
system and thus of its adjoint ([T|). The Krein-Rutman theorem gives the n^'s, but the solution is not 
explicitly known except for / = 1, a situation where the motor effect cannot be achieved. The stability 
of this problem is also related to the notion of relative entropy [El [211 [221 [20] . 



The typical results obtained about biomotors in [71 [16] are that, for small diffusion a, under some 
precise asymmetry assumptions on the potentials, the solutions tend to concentrate, as o" ^ 0, as 
Dirac masses at either x = or x = 1. In the sequel such a behavior will be called motor effect. 

Our results (i) provide an alternative proof of this motor effect, and (ii) allow for more general 
assumptions like, for instance, various scalings on the coefficients vij. While [3 [16] transform the 
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system ([T]) into an ordinary differential equation and analyze directly its solution, here we use a 
direct PDE argument based on the phase functions Ri = —a In rii that satisfy (in the viscosity sense, 
[21 El El [H] ) a Hamilton- Jacobi solution. This is reminiscent to the method used for front propagation 
( \13\ S]). We recall that the appearance of Dirac concentrations in a different area of biology (trait 
selection in evolution theory) relies also on the phase function and the viscosity solutions to Hamilton- 
Jacobi equations, [TOl [5], 

In Section [2] we obtain new and more precise versions of the results of [7J by analyzing the asymp- 
totics/rates as o" ^ 0. In Section [3] we present new results for large transition coefficients, while in 
Section H] we consider coefficients that may vanish. 

2 Bounded non-vanishing transition coefficients 

We begin with the assumptions on the transition rates and potentials. As far as the former are 
concerned we assume that 

there exists k > such that fjj > A; > for all i ^ j. (6) 

As far as the potentials are concerned we assume that, for all i = 1,2, 

i^iGC^'\0,l), (7) 
there exists a finite collection of intervals {Jk)i<k<M such that mvnip^yQ in IJjfc, (8) 

and 

maxV''>0 in [0,11. (9) 

l<i<I 

Notice that these assumptions are satisfied by periodic potentials with period 1/M. 
Figure 1: Motor effect exhibited by the parabolic system dlj with two asymmetric potentials. Left: 

THE potentials Ipi, tp2- RiGHT: THE PHASE FUNCTIONS Rf = — crln(ni), R2 = ~CTln(n5). AS ANNOUNCED IN 

Theorem 12. 11 we have Ri x R2 and are nondecreasing. This means that the densities are concentrated 

AS DiRAC MASSES AT a; = 0. HERE WE HAVE USED a = 10"''. SeE FIGURE [2] FOR ANOTHER BEHAVIOR. 

Our first result is a new and more precise version of the result in [7]. It yields that the system ([T]) 
exhibits a motor effect for a small enough and molecules are necessarily located at x = 0. This effect 
is explained by a precise asymptotic result in the limit cr — > 0. 

To emphasize the dependence on the diffusion a, in what follows we denote, for all i = 1, by 

the solution of ([T|). Moreover, instead of we use the normalization 

l<i<I 

We have: 

Theorem 2.1 Assume that (0), ^ and (C^j hold. Then, for all i = 

RI = -a\nn1 R in C(0, 1), i?(0) = and R' = miii^('(/'-)+ ■ 
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In physical terms, R can be seen as an effective potential for the system. To state the next result, 
we recall that throughout the paper we denote by 5o the usual (5-function at the origin. 
We have: 

Corollary 2.2 Assume, in addition to ^ and that min ■0^(0) > 0, and normalize 

nf by ^ instead of fTOj) . There exist {pi)i<i<i such that 

Pi^o, Pi > 0, and ^ Pi = I- 

l<i<I 



There are several possible extensions of Theorem I2.1[ Here we state one which, to the best of our 
knowledge, is not covered by any of the existing results. 

To formulate it, we need to introduce the following assumption on the potentials {ipi)i<i<i which 
replaces ([9]) and allows to consider more general settings. It is: 

the set {x G [0, 1] : max tp'iix) < 0} is a union of finitely many intervals (-f^Oi<KM') 

!<«</ 

< 

and (UJ^y n {UKiY is either a finite union of intervals or isolated points. 

(11) 

We have: 

Theorem 2.3 Assume (Qj, {TiP and figj. Then 

R1 = -C7ln< R in C(0, 1), i?(0) = and 

^min^(''Ai)+ in ^ Jk, 



R'=l 



max ibi in U Ki, 
i<i<l 



^ m lnt{{UJkrn{UKiY). 



As a consequence we have: 



Corollary 2.4 In addition to {3), {HP o-nd assume that we have the same number 

of intervals and Ki in ^ and Ml]) respectively, that is the left endpoint of Ji and, finally, that, 
for all k = l, ...,M, 



max ■ip'^{y)dy\ < / min ■ip'^{y)dy. 
Then, for all i = 1, there exist {pi)i<i<i such that 

rii Pi6o, Pi > 0, and ^ Pi = I- 

l<i<I 
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Other possible extensions concern coefficients that may vanish somewhere and/or be unbounded. 
The former case is studied in Section HI As far as the Vij being unbounded, it will be clear from the 
proof of Theorem 12 .H that the coefficients can depend on a as long as, for a ^ and all i, j = 1, I, 
there exists a > such that 







and 



oo. 



Going further in this direction leads to a different limits for — In nf that we study in the next Section. 



We continue next with the proof of Theorem 12. li The modifications needed to prove Theorem 
are indicated at the end of this section where we also discuss the proofs of the Corollaries. 



Proof of Theorem 12.11 A direct computation shows that the -Rf s satisfy, for va = vi,, 
system 



V^-', the 



dx 



avii in (0, 1) 



dx 



i^'i in {0, 1}. 

Adding the equations of ([T]) and using ([3]) yield the conservation law 

-Si E <i - ^1 E = o. 

l<i<I l<i<I 

which together with the boundary condition gives 

d 



l<i<I 



Setting 



we have 



l<i<I 



dx 

and, as a consequence, the total flux estimate 



1 

<i<i 



mm < — — < max ip^. 



l<i<I 



dx 



l<i<I 



fl2) 



(13) 



(14) 



The normalization (jlOp of the n^'s implies that S"^(0) = 0. As a result, there exists a S £ C ' (0, 1) 
such that, after extracting a subsequence, 



S, S{0) = 0, and 



dS 



(15) 



min ih'^ < -— < max ib!: in \0,1\. 
i<i<i dx ~ i<i<i 



Next we obtain bounds on the -Rf s, which are independent of a, and imply their convergence as 
cr — > 0. This is the topic of the next Lemma which we prove after the end of the ongoing proof. 
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Lemma 2.5 For each i = 1,...,/ there exists a positive constant Ci = Ci{ip^, ai^u, aip'-)) such that 

. dm . 



Moreover, for all i = 1, 



Ri\ + \^\<a in [0,1], 



R^~R = S, mC([0,l]). 



We obtain next the Hamilton-Jacobi satisfied by the limit R = S. The claim is that the limit is a 
viscosity solution (see, for instance, [3l|9]) of 



OR 



dx 



+ max M,'— ] = in (0,1). 

i<i<i ox 



(16) 



We do not state the boundary condition because we do not use them. It can, however, be proved 
that R satisfies 

dR , dR 

— — < max at x = and — — > mm ih, at x = 1. 

dx i<i<i dx l<i<I 

We begin with the subsolution property. Letting cr ^ in the inequality 



dx"^ 



-a- 



dRI 



dx 



^ives, for all i = 1, /, 



dR 



dx 



,,dR 



To prove that i? is a supersolution of (jl6p we observe that function R" = min^RI satisfies the 



inequality 



d^ 

dx"^ 



-a- 



dR° 



dx 



dR"^ 



+ max [-i^'i{x)-^] + 0- ^ T^ij > o"min(^'jj). 

«J=1 



Letting again a — > 0, we find that R = S = lim R'^ satisfies 



dR 



dx 



. 1 1 / \ dR, 
+ max \—Wi{x)^—\ > 0. 
i<i<i ' dx 



We obtain now the formula for R. To this end, observe first that, since lim R^ = R = S, letting 



(T ^ in ([M]) yields 



dR 

min ib'i < — — < max ■0'. (17) 
i<i</ * ~ dx ~ i<i<i 



Next we show that, in the viscosity sense. 
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Indeed for a test function let xq € (0,1) be the maximum of i? — i.e., {R — $)(xo) = 
max (R — and assume that 

0<a:<l 

$'(xo) < 0. 

Applying the viscosity subsolution criterion to (I17|) . then implies that 

$'(xo) - maxV'i(xo) > 0. 

i 

This, however, contradicts the inequality 

max ^Uxq) > 
l<i<7 

that follows from the assumption ([8]). 
Combining (fTTl) and (fTSj) we get 

min < < ^ax ■0,'. (19) 

Finally, given a test function let xq G (0, 1) be such that {R — ^){xq) = max {R — <I>) and 

0<a;<l 

assume that 

$'(xo) > 0. 

Again by the viscosity criterion we must have 



$'(xo) — min ^pi{xo) < 0, 

!<«</ 



and, hence, in the viscosity sense. 



dR , ,,, „ dR , ^ 

— < min iPl)+ if 1^ > 0. 20 

dx i<i<i dx 

This concludes the proof of the formula in the claim. □ 
We return now to the 

Proof of Lemma 3.1 For the Lipschitz estimate, observe that, at any extremum point xq of -q^, 
we have = 0. Evaluating the equation at xq, we get 



dRI ^ . ^,dR1 



As a consequence, at xq we have 



dx 



dRI 



< max + V crm- 

0<x<l 



dx 

To identify the limit of mini<j</ Rj notice that the inequality 



n-i S > n,- s J maxn^ 
i<i</ ^ 
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gives 

-C7ln/+ mill m <S" < R1, 



and thus 



S" < min Rf. 

i<i<l 



Consequently, we have the uniform convergence 

min Rf rr* S. 

l<i<I * '^^0 

To prove the claim about the limit of the Rf we observe that summing over i the equations of (fT2 
yields 

i,i=l ij=l !<«</ 

Integrating in x and using the gradient estimates, we find that 



RDl < Co. 



Together with the uniform gradient estimate on R^ and the uniform bound on mini<j</i?J, we 
deduce that 

□ 

We continue with the 

Proof of Corollary 12.21 The normalization (jll) amounts to adding a constant to the Ri. The 
exponential behavior of n?", with an increasing R^ (from Theorem 12. ip . yields that the nf s converge, 
as (7 ^ 0, to uniformly on intervals [e, 1] with e > 0. Moreover, i2(0) = 0. The result follows with 
Pi > 0. If yOj = for some i = 1, /, then, letting o" ^ in ([1]), gives, in the sense of distributions, 
that 

But then all the pj must vanish, which is impossible with the normalization of unit mass. □ 
We present now a brief sketch of the proof of Theorem 12.31 Since it follows along the lines of the 
proof of Theorem l2.lt here we only point out the differences. 
We have: 

Proof of Theorem 12.31 The Lipschitz estimates, the passage in the limit and the identification of 
the limiting Hamilton- Jacobi equation in the Theorem 12.11 did not depend on the assumption ([9]) , 
hence, they hold true also on the case at hand. The final arguments of the proof of Theorem 12.31 also 
identify the limit on the set (UX;)'^. On the set U-fC/ we already know from (jl7p that R! is less than 
the claimed value, and thus it is negative. We conclude the equality by using the Hamilton-Jacobi 
equation. Indeed in this situation we know that 

dR dR 
max \—^pi^:^\ = —-^ max w,-. 
i<i<i dx dx i<i<i 
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□ 

We conclude the section with the proof Corollary 2.4, which is simply a variant of the one for 
Corollary 2.2. We have: 

Proof of Corollary 2.4 The assumption on UJ asserts that R is increasing on UJ. Then it may 
decrease but, for x > 0, R{x) > R{0). With the unit mass normalization, this means that R{0) = 
as before and the convergence result holds as before. □ 



3 Large transition coefficients 



Figure 2: Motor effect exhibited by the parabolic system ([2T|) with large transition coefficients. The 

FIGURE DEPICTS THE PHASE FUNCTIONS i?i , i?f . AS ANNOUNCED IN THEOREM 13.11 WE HAVE Ri f» _Rf AND CAN 
DECREASE SLIGHTLY. HeRE WE HAVE USED a = 5 10^^. 

In this section we consider transition coefficients normalized by 1/cr. For the sake of simplicity we 
take 1 = 2. This allows for explicit formulae. The equations for larger systems, i.e., I > 3, are more 
abstract. The system ([T]) is replaced by 

' - <) + iz^m? = ii.2n^ in (0, 1), 

' -^-^ri<^2-U^^2nl) + lv2n- = lv^n'l in (0,1), (21) 

^ + VV'i nf = in {0,1} for i = l,2. 

As before we assume that 



< > in [0,1] for i = 1,2. 



(22) 



The result is: 

Theorem 3.1 Assume (0), (0), (Q), (0) and (0) and consider the solution {n'l^n'^) to \21\) normalized 
by (0) + (0) = 1. Then, as a ^ and i = 1, 2, 



RI = -a\nn1 R in C(0, 1), i?(0) = 0, anc 



R' > < 



min lb': on U J/, 
i<i<r' 



^ -Vk on {UJiY 



The corollary below follows from Theorem 13.11 in a way similar to the analogous corollaries in the 
previous section. Hence, we leave the details to the reader. 

Corollary 3.2 In addition to the assumptions of Theorem \3.1{ suppose that G Ji, the potentials are 
small enough so that 

Vk\K\< / mm v)Av)av, 

J sup J l<i<2 
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and {711,112) is normalized by There exist pi,P2 > such that /?i + /02 = 1 CLn-d, as o" — > and for 
i = l,2, 

We present next a sketch of the proof of Theorem l3.1l as most of the details follow as in the previous 
theorems. 

Proof of Theorem l3.1I The total flux and Lipschitz estimates follow as before. The main new point 
is the limiting Hamilton-Jacobi equation which is more complex. We formulate this as a separate 
lemma below. Its proof is based on the use of perturbed test functions. We refer to [4] for the rigorous 
argument in a more general setting. 



Lemma 3.3 The uniform in [0, 1] limit R, as a ^ 0, of the Rf satisfies the Hamilton-Jacobi equation 

dR 

■ dx 



OR 

H{—,x)=0, %n (0,1) (23) 



with 

Hip, x) = i [f3i + f32 + V + 02? -mP2- 1^11^2] , (24) 

where, for i = 1,2, 

The formula for R' follows from the above Lemma by analyzing the solutions to the Hamilton- 
Jacobi equation as before. On the set UJi the answer follows from the bounds (|14p . On the set {Ujy 
the argument is more elaborate. Using that R' is a subsolution, we get 

P1P2 - viV2 > and Pi+p2< 0. 
Therefore both Pi and P2 are nonpositive and thus 

{R'f - i)[R! -Vi<^. 

On the other hand we know that on {^JiY one of the potentials - for definiteness say ipi - satisfies 
■0^ > 0, hence, always in (UJ/)'^, 



ii'>-W(x)-^(V'i)2 + 4z.i]> V^. 
The inequalities for R! are now proved. □ 

4 Vanishing transition coefficients 

We focus here to the case where the transition coefficients {i^ij)i<ij<i may vanish at either some points 
or, in fact, on large sets. In this situation, we assume that 

for each j = 1, ■0^ < on a finite collection of intervals {Kj')i<a<Aj and 

for all j = 1, / and a = 1, .■.,Aj, there exists i £ {= 1, /} such that (25) 
> on Kj^, and, in a left neighborhood of the right endpoint of Kj', Uij > 0. 
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To go for weaker assumptions would face the completely decoupled case (when v vanishes) and the 
motor effect does not occur. 



We have: 

Theorem 4.1 Assume (0), (0^, (E^i o^'^d normalize the solution (nf)i<j</ to [21\) by [TU) . For 
i = I, /, let R1 = — Innf . Then, as a ^ 0, 

either Ri mC(0, 1), or R^ ^^^> oo uniformly in [0,1]. 

Moreover, the function R = min^ Ri satisfies 

R(0) =0, R' >0 and R' = min ^h' on U J,. 

We also have: 

Corollary 4.2 In addition to the assumptions of Theorem \^.l\ suppose that G Ji. For i = 1,,,/, 
there exist pi > such that ^ki^ci Pi — and, as a ^Q, 

The direct conclusion of Theorem 14.11 is simply that 

l<i<I 

The corollary follows from the fact that, for all i = 1, > 0. We do not know whether in 

this context each pi is positive. To get this, we need to assume something more like, for example, 
fij(O) > for ah i,j = 1, /. 

We conclude with a brief sketch of the 
Proof of Theorem 14. li The total flux and Lipschitz estimates follow as before. A careful look 
at the proof of the convergence part of Theorem 12.11 shows that either the s blow up or they are 
uniformly bounded and, hence, converge uniformly in (0, 1) to a subsolution of 

\R',\^-4R',<0. 

It then follows that 

i?. > on (UaK^'y and R- < ^-(x) on U„ iCf . 
The final step is to prove that 

R{x) = min Ri{x) in L{x) = {i, ipi^x) > 0}. 

ieL{x) 

This follows as before. We leave the details to the reader. □ 
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